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ABSTRACT 

Although many long non-coding RNAs (IncRNAs) 
have been discovered, their function and their asso- 
ciation with RNAi factors in the nucleus have 
remained obscure. Here, we identify RNA tran- 
scripts that overlap the cyclooxygenase-2 (COX-2) 
promoter and contain two adjacent binding sites 
for an endogenous miRNA, miR-589. We find that 
miR-589 binds the promoter RNA and activates 
COX-2 transcription. In addition to miR-589, fully 
complementary duplex RNAs that target the COX-2 
promoter transcript activate COX-2 transcription. 
Activation by small RNA requires RNAi factors 
argonaute-2 (AG02) and GW182, but does not 
require AG02-mediated cleavage of the promoter 
RNA. Instead, the promoter RNA functions as a 
scaffold. Binding of AG02 protein/small RNA 
complexes to the promoter RNA triggers gene acti- 
vation. Gene looping allows interactions between 
the promoters of COX-2 and phospholipase A2 
(PLA2G4A), an adjacent pro-inflammatory pathway 
gene that produces arachidonic acid, the substrate 
for COX-2 protein. miR-589 and fully complementary 
small RNAs regulate both COX-2 and PLA2G4A gene 
expression, revealing an unexpected connection 
between key steps of the eicosanoid signaling 
pathway. The work demonstrates the potential for 
RNA to coordinate locus-dependent assembly of 
related genes to form functional operons through 
c/s-looping. 



INTRODUCTION 

Many loci encode long non-coding RNAs (IncRNAs) that 
overlap the 5' and 3' termini of genes (1). Although the 
functional significance of many annotated nuclear 
IncRNAs remains under debate (2), recent reports have 
implicated them in the control of transcription (3-5). 
Understanding how IncRNAs may regulate gene tran- 
scription would provide insights into a mechanism that 
has the potential to supplement control by protein tran- 
scription factors. In addition, IncRNAs that regulate ex- 
pression of key genes would be new targets for therapeutic 
intervention. 

Well-studied examples of IncRNAs in mammalian cells 
include XIST and TSIX, RNAs that control X-chromo- 
some inactivation (6), and HOTAIR, a IncRNA that con- 
tributes to silencing at the HOX locus (7). These RNAs 
appear to alter gene expression by functioning as scaffolds 
for histone modification complexes. Published models, 
however, cannot explain: (i) the basis for recognition at 
a specific gene, (ii) how small RNAs in concert with RNAi 
factors affect gene transcription, and (iii) how an IncRNA 
mediates sequence-specific control of gene expression. 

LncRNAs functioning in trans induce effects at spatially 
distant loci. In contrast, action in cis is restricted to genes 
that are spatially close to where the IncRNA is 
transcribed. Localization of c/s-acting IncRNAs to the 
chromatin near their sites of transcription may allow for 
the regulation of nearby genes. Cis interactions could, 
therefore, mediate a new layer of gene-specific regulation 
by complementary small RNAs and RNA-binding 
proteins to alter transcription, splicing, histone modifica- 
tion or local chromatin structure. 
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To test whether c/s-acting RNAs can alter gene tran- 
scription, we focused on the genomic locus encoding 
cyclooxygenase-2 (COX-2) and phospholipase A2 group 
IVA (PLA2G4A) (Figure 1A). COX-2 protein is a critical 
regulator of inflammation in normal physiology and 
disease. COX-2 catalyzes the conversion of arachidonic 
acid to intermediates that are further metabolized to pros- 
taglandins and other eicosanoids (8). These eicosanoids 
mediate numerous biological processes including inflam- 
mation (9), development (10), reproduction (11), cellular 
immunity (12), cancer (13) and energy homeostasis (14). 
The COX-2 promoter contains a TATA-box and 
numerous binding sites for transcription factors (15). 
Like many pro-inflammatory genes, transcription factors 
are recruited to these sites in different combinations de- 
pending on which signaling pathways are activated (16). 

PLA2G4A is a cytosolic phospholipase that catalyzes 
the hydrolysis of membrane glycerophospholipids to 
release arachidonic acid, the substrate for the COX-2 
enzyme (17). The gene encoding PLA2G4A is located on 
chromosome lq25, adjacent to the COX-2 locus. The 
PLA2G4A and COX-2 promoters are organized head to 
head but are separated by 149 000 bases. The PLA2G4A 
promoter lacks SP1 sites and stretches of CpG islands that 
are characteristic of TATA-less promoters (18), suggesting 
that its expression may be influenced by other regulatory 
factors and pathways. Although the functions of COX-2 
and PLA2G4A proteins are intimately linked, it is not 
known how their expression may be coordinated and 
how co-regulation of these pro-inflammatory genes 
would impact eicosanoid signaling. 

The importance of COX-2 in normal physiology, 
disease and therapeutic development makes it a high-pri- 
ority target for understanding gene regulation. Here, we 
address the role of RNA-directed interactions in 
regulating COX-2 and PLA2G4A transcription. Unlike 
protein transcription factors that are general regulators 
of many genes, the selectivity inherent in RNA recognition 
permits sequence-selective control of gene targets. We dis- 
covered a network of RNAs overlapping the COX-2 
promoter leading us to investigate their impact on the 
expression of COX-2 and PLA2G4A. Both genes are 
activated by endogenous miR-589 that targets sites 
within a COX-2 promoter RNA, demonstrating the po- 
tential for RNA to coordinate transcription of genes 
within the eicosanoid signaling pathway. 

MATERIALS AND METHODS 

mRNA-seq library generation and sequencing 

A549 cells were first harvested, and pure nuclei were 
isolated. The nuclei were then dissolved in TRIzol, and 
total nuclear RNA was isolated subsequently following a 
standard TRIzol RNA preparation protocol (Sigma). 
Total nuclear RNA was treated with DNase I 
(Worthington) for 20min at 37°C. 

Illumina TruSeq mRNA-seq libraries were constructed 
according to manufacturer's instructions (Illumina, Inc., 
San Diego, CA, USA) by starting with 1 jig of total RNA. 
To remove the most abundant duplex DNA species 



(representing more highly expressed transcripts), the 
mRNA-seq library was subjected to duplex-specific 
nuclease (DSN) treatment. In all, 280 ng of the TruSeq 
mRNA-seq library was treated with DSN according to 
manufacturer's instructions (Evrogen/Axxora, LLC, San 
Diego, CA, USA) using a 2-h annealing period before 
DSN addition. Following DSN treatment, the library 
was subjected to an additional 14 cycles of PCR and 
post-PCR clean up using Agencourt AMPure XP beads 
(Beckman Coulter, Brea, CA, USA). 

RNA-seq libraries were sequenced using the Illumina 
Hiseq 2000 as per manufacturer's instructions for 
paired-end 2x 100. In all, 108 million raw paired reads 
(100 nt) were obtained and then aligned to human 
genome hgl9 using TopHat 1.4.1 (and Bowtie 0.12.8) by 
allowing up to two mismatches to the reference sequence 
in a paired-end mapping mode. Approximately, 71% of 
the raw reads were uniquely aligned (-g 1). The aligned 
reads were visualized in Integrative Genomics Viewer 
program. 

Small RNA-seq library generation and sequencing 

A549 cells were harvested, and pure nuclei were isolated. 
The nuclei were then dissolved in TRIzol, and total 
nuclear RNA was isolated subsequently using Qiagen's 
miRNeasy Mini Kit. Total nuclear RNA was treated 
with DNase I (Worthington) for 20min at 37°C. 

The small RNA-seq library was made following the 
standard Illumina small Truseq preparation protocol. 
RNA-seq libraries were sequenced using the Illumina 
Hiseq 2000 as per manufacturer's instructions for single- 
end 1 x 50. Roughly, 50 million raw paired reads 
(averaged over two replicates) were obtained. Each raw 
read was trimmed by removing the adaptor sequence. 
The trimmed reads shorter than 15nt were removed. 
The filtered trimmed reads were aligned to human 
genome hgl8 using Bowtie 0.12.7 by allowing up to two 
mismatches to the reference sequence. Up to 10 different 
alignments per read was permitted. For reads that were 
aligned to multiple positions in the reference genome, the 
single aligned read with the fewest mismatches was 
selected using a Perl script. If there were still more than 
one aligned reads in this case, the original read would be 
disregarded. Approximately, 72% of the raw reads were 
successfully aligned. The aligned reads were again mapped 
to UCSC miRNA database and miRBase (mature 
miRNA) to search for possible miRNA hits. 

Rapid amplification of cDNA ends 

Rapid amplification of cDNA ends (RACE) was per- 
formed using the Gene Racer Kit (Invitrogen). cDNA 
was prepared from total RNAs of A549 cells through 
reverse-transcription (RT) reaction using oligo dT or 
random primers. The 5' or 3' end of cDNA was amplified 
through two PCRs using Platinum Taq DNA polymerase 
(Invitrogen) and multiple primer sets (Supplementary 
Table SI). The thermal cycling condition for the first 
PCR was 94° C for 2min, followed by 5 cycles of 94° C 
for 30 s and 72° C for 1 min, 5 cycles of 94° C for 30 s and 
70°C for 1 min, and 25 cycles of 94°C for 30 s, 66°C for 
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Figure 1. Transcription at the COX-2 locus is bidirectional and coordinated. (A) Eicosanoid biosynthetic pathway. (B) Schematic of transcripts at 
the COX-2 promoter identified by 5 r and 3 r RACE. (C) Identification and quantification of transcriptional domains at the COX-2 promoter by 
strand-specific RT-qPCR. n = 3. (D) Expression of COX-2 mRNA or promoter RNA in cancer cell lines, n = 3. (E) Expression of COX-2 mRNA or 
promoter RNA in human tissues. (F) Expression of COX-2 mRNA or promoter RNA in myometrial tissue from patients. (G, H) Effect of adding 
pro-inflammatory agents on levels of COX-2 (G) protein, (H) mRNA and promoter RNA. Vehicle is solvent only control. Cells were treated with 
epidermal growth factor (lOng/ml), PMA (50ng/ml), IL1J3 (lOng/ml) or TNFa (lOOng/ml) for 24 h. n = 3. Error bars are SD. Experiments use A549 
cells unless otherwise noted. Westerns are representative of multiple experiments. See also Supplementary Figure SI. 
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30 s and 68° C for 1 min. The condition for the second 
nested PCR was 94°C for 2 min, followed by 25 cycles 
of 94°C for 30 s, 65°C for 30 s and 68°C for lmin. 
PCR products were analyzed on 1% agarose gel 
(Supplementary Figure SIB). Major PCR products on 
gels were cloned and sequenced to determine transcription 
start sites or 3' ends of transcripts (Supplementary 
Figure SIC). 

Quantitative PCR/strand-specific RT-qPCR 

Total RNA from human cells and tissues were isolated 
using Tri-Reagent (Sigma). Two micrograms of total 
RNA were treated with DNase I to remove genomic 
DNA. Treated RNAs were reverse transcribed using the 
High Capacity cDNA Reverse transcription kit (Applied 
Biosystems). Random or strand-specific primer (exon 3 R 
or —740; 0.5 uM) was used in the RT reaction. 
Quantitative PCR (qPCR) was performed using iTaq 
SYBR Green Supermix (Biorad) with 50 ng of cDNA as 
template. Data were normalized relative to measured 
GAPDH levels. Primers used in qPCR are listed in 
Supplementary Table SI. 

Cellular delivery of duplex RNAs, single-stranded gapmers 
and locked nucleic acid miR inhibitors 

Lipofectamine RNAiMAX (Invitrogen) was used in 
reverse transfection experiments (i.e. lipid/RNA mixtures 
added to a culture well before cells) to deliver duplex 
RNAs, single-stranded gapmers or locked nucleic acid 
(LNA) anti-miRs into A549 cells as described (19). Cells 
transfected with duplex RNAs were harvested 3 days after 
transfection for qPCR and 4-5 days after transfection for 
western blot. In double transfection experiments with 
duplex RNAs, the second transfection was performed 2 
days after the first transfection, and cells were harvested 
on day 5. Sequences of duplex RNAs and gapmers are 
listed in Supplementary Tables S2 and S3. Anti-miRs 
complementary to miR- 5 89 and -452 were obtained from 
Exiqon. Transfection of single-stranded gapmers and 
LNA anti-miRs required a double transfection. Reverse 
transfections were performed in 6-well plates using 125k 
A549 cells/well. Three days after the first transfection, cells 
were counted and reverse transfected a second time. Cells 
were harvested 3 days later for western blot and qPCR. 
Concentrations for each experiment are provided in the 
figure legends. Transfection of T47D required 200k cells/ 
well. 

Western analysis 

Cells were harvested with trypsin-EDTA solution 
(Invitrogen) and lysed. The protein concentration in 
each sample was quantified by BCA assay (Thermo 
Scientific). SDS-PAGE was performed using 7.5% Tris- 
HC1, 7.5% TGX or 4-20% TGX gels (Biorad). Gels were 
run at 100 V for 70 min. After gel electrophoresis, proteins 
were transferred to nitrocellulose membrane (Hybond-C 
Extra; GE Healthcare Life Sciences) at 100 V for 90 min. 
Membranes were incubated with specific primary 
antibodies at the following dilution ratios: anti- 
COX-2 antibody (1:2500; 160112; Cayman Chemical), 



anti-cPLA2 (PLA2G4A) (1:2000; 5479S; Cell Signaling), 
anti-P-Actin (1:20 000; A5441; Sigma), anti-AG02 
(1:1000; ab57113; Abeam) and anti-GW182 (1:4000; 
A302-329A; Bethyl Labs). Horseradish peroxidase- 
conjugated anti-mouse (715-035-150) or anti-rabbit (711- 
035-152) secondary antibody (1:4000-1:20 000; Jackson 
Immunolabs) was used for visualizing proteins using 
SuperSignal West Pico Chemiluminescent Substrate 
(Thermo Scientific). 

Melting temperature ( J m ) determination 

Annealed duplex RNAs were diluted to 1.5 uM in 0.1 M 
NaH 2 P0 4 /Na 2 HP0 4 buffer (pH 7.4) (final volume 400 ul) 
and transferred to quartz cuvettes. The solutions were 
topped with mineral oil (300 jal) and heated from 15 to 
95° C at l°C/min while monitoring changes in A260 
using a Cary 100 Bio spectrophotometer. Samples were 
then re-annealed in the spectrophotometer, and the experi- 
ment was repeated several times. T m values were 
calculated as the first derivatives of the plot of A260 
versus temperature, and the results from three to four 
experiments were averaged. 

Treatment of untransfected/transfected cells with 
physiologic agents 

Two days after plating, cells were treated with epidermal 
growth factor (EGF; lOng/ml), phorbol 12-myristate 13- 
acetate (PMA; 50ng/ml), interleukinl-beta (IL1(3; 10 ng/ 
ml), tumor necrosis factor-alpha (TNFoc; 100 ng/ml), dexa- 
methasone (Dex; 1 uM) or vehicle (water or DMSO) in 
media containing 2% charcoal-stripped FBS (Atlanta 
Biologicals) for 24 h then harvested for qPCR or western 
blot. For combination treatment with RNA 12 (25 nM) or 
RNA12ftc (25 nM), compounds were added to cells 2 days 
after transfection in media containing 2% charcoal- 
stripped FBS and harvested 24 h later for RNA. qPCR 
was used to evaluate levels of mRNA and IncRNA. 

Chromatin immunoprecipitation/RNA 
immunoprecipitation 

Chromatin immunoprecipitation (ChIP) and RNA 
immunoprecipitation (RIP) assays were performed as pre- 
viously described (20,21). A549 cells (2 000 000 cells in 
15-cm dish) were reverse-transfected with RNA 12, 
RNA12^c or mismatch oligomer (25 nM) on Day 0. 
Untreated cells were also prepared on Day 0. The cells 
were treated with either vehicle or 10 ng/ml IL1(3 for 1 h 
on day 3. For most ChIP and RIP cells were crosslinked 
with 1% formaldehyde for 10 min, whereas no crosslink 
protocol (20) was used for AG02-RIP, and 0.1% formal- 
dehyde was used for GW182-RIP. Cells were recovered by 
scraping, and nuclei were isolated using hypotonic lysis 
(HL) buffer [4ml; 10 mM Tris-HCl (pH 7.5), 10 mM 
NaCl, 3mM MgCl 2 , 0.5% (v/v) NP-40]. Nuclei were 
lysed in lysis buffer [1ml; 1% SDS, 10 mM EDTA, 
50 mM Tris-HCl (pH 8.1), lxRoche protease inhibitors 
cocktail, RNasin Plus RNase inhibitor (Promega)] and 
sonicated (2 pulse, 20% power setting, 20 s). 

Nuclear lysate (40-250 ul) was incubated overnight 
with specific antibodies (2-5 ug) in immunoprecipitation 
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buffer. After antibody-protein-DNA/RNA complex was 
recovered with 40-50 ul of protein G plus/protein A 
agarose beads (Calbiochem), the beads were washed 
with buffers. The complex was eluted twice with 250 ul 
of elution buffer (1% SDS, 0.1 M NaHC0 3 ). Cross- 
linking was reversed by adding NaCl to a final 
concentration of 200 mM and heating at 65°C for at 
least 2h. Protein was digested by incubating with protein- 
ase K (20 ug; invitrogen) at 42°C for 50min, followed by 
phenol-chloroform extraction and ethanol precipitation. 
For ChIP, samples were dissolved in 200 ul of nuclease 
free water and analyzed by qPCR. For RIP, samples 
were dissolved in 17.6 ul of nuclease free water and 
treated with DNase I at 25°C for 20min. RT reactions 
(20 ul) were performed for input and +RT samples, then 
the samples were analyzed by qPCR. PCR products were 
analyzed by 2.5% agarose gel electrophoresis, and gels 
were stained with ethidium bromide. 

Antibodies used in ChlP/RIP were as follows: anti- 
RNAP2 (05-623, Millipore), anti-H3K4me3 (ab8580, 
Abeam), anti-H3K27me3 (07-449, Millipore), anti-H4Ac 
(06-598, Millipore), anti-NFucB p65 (17-10 060, Millipore), 
anti-NFKB p50 (06-886, Millipore), anti-CREBl (17-600, 
Millipore), anti-AG02 (015-22 031, Wako), anti-WDR5 
(A302-429A, Bethyl Labs) and anti-GW182 (A302-329A, 
Bethyl Labs). Normal mouse immunoglobulin G (IgG) 
(12-371, Millipore) or normal rabbit IgG (12-370, 
Millipore) was used as a negative control. The sequences 
of qPCR primers used in ChIP or RIP are listed in 
Supplementary Table S 1 . 

Cell fractionation and extract preparation 

A549 cells were transfected with 25 nM RNA12, 
RNA12ftc or mismatched control duplex (MM) using 
Lipofectamine RNAiMAX and harvested 24 h later. 
Cells were fractionated into cytoplasmic and nuclear 
extracts. Cells were first resuspended in HL buffer 
[10 mM Tris (pH 7.4), 10 mM NaCl, 3mM MgCl 2 , 0.3% 
NP-40, 1 % protease inhibitor cocktail] then incubated on 
ice for 15min and the soluble fraction kept as cytoplasmic 
extract. To the cytoplasmic extract, KC1 and glycerol were 
added to a final of 140 mM and 10%, respectively. Nuclei 
in the pellet were washed 3x with HL buffer then resus- 
pended in either nuclear lysis buffer [20 mM Tris (pH 7.4), 
0.15M KC1, 3mM MgCl 2 , 0.5% NP-40, 10% glycerol, 
1% protease inhibitor cocktail]. Nuclei resuspended in 
nuclear lysis buffer were sonicated on ice at 20% power 
3x for 20 s. The soluble fraction was kept as nuclear 
extract. All extracts were stored at — 80°C for later use. 

Samples from whole cell, cytoplasmic and nuclear frac- 
tions were analyzed by western blot to show cellular lo- 
calizations of AG02 and GW182. Western analysis was 
performed using anti-AG02 (1:1000; ab57113; abeam), 
anti-GW182 (1:4000; A302-329A; Bethyl Labs), anti- 
tubulin (1:6000; T5201; Sigma), anti-calreticulin (1:1000; 
2891S; Cell Signaling) or anti-histone H3 (1:10 000; 2650S; 
Cell Signaling) primary antibody and horseradish perox- 
idase-conjugated anti-mouse IgG (715-035-150) or 
anti-rabbit IgG (711-035-152) secondary antibody 



(1:5000-1:10 000; Jackson Immunolabs) as described 
earlier in the text. 

Co-immunoprecipitations 

Co-immunoprecipitation experiments were performed by 
mixing 30 ul of Protein G Plus/Protein A resin, 1.5 jig of 
antibody (anti-AG02, ab57113, Abeam; anti-GW182, 
Bethyl Laboratories) and nuclear extract (0.5-1 mg of 
total protein) and incubating at 4°C for 4h. Resin was 
washed 4x with IP wash buffer [20 mM Tris (pH 7.5), 
0.4 M NaCl, 2mM MgCl 2 , 0.05% NP-40, 0.025% SDS] 
and co-purified proteins eluted by boiling resin in 25 ul 1 x 
SDS loading buffer. Elution was resolved by SDS-PAGE, 
and proteins detected by western blot with anti-AG02 or 
anti-GW182 antibodies. 

AG02 cleavage assay 

AG02 protein was immunoprecipitated from each extract 
using 1.5 jig of anti-AG02 antibody (Abeam, ab57113) 
and Protein G Plus/Protein A agarose. Target RNA sub- 
strate RNAllSub was 5' radiolabeled and gel-purified. 
Labeled RNAllSub (lOOOOOcpms, M).2pmols) was 
added to the resin in 1 x RNAi buffer [50 mM Tris (pH 
7.4), 2mM MgCl 2 , 0.5 mM DTT, 100 mM KC1, 50 mM 
NaCl, ImM EDTA] supplemented with 0.5ug/ul yeast 
tRNA and 2.5 mM ribonucleoside triphosphates. 
Cleavage reactions were incubated at 30°C for 1.5 h, and 
RNA was collected by phenol extraction. Extracted RNA 
was precipitated with 9 volumes of 2% LiC10 4 in acetone, 
washed with acetone, resuspended in 90% formamide, 
lx TBE buffer and resolved on a 14% denaturing poly- 
acrylamide (7M urea, lx TBE). Gels were dried and 
exposed to phosphorimager screens overnight to visualize 
radioactivity. KNAllSub consists of a single siRNA 
target site flanked by firefly luciferase gene sequence and 
capped with two deoxynucleosides on each end. 
RNA \2Sub sequence (siRNA target site capitalized): 
gaacaauugauuuuacagacGACAAUUGGUCGCUAACC 
Gucacguacgcggaauacutc . 

Chromosome conformation capture (3C) assay 

A549 cells (2 million cells per 15-cm dish) were reverse- 
transfected with RNA 12 or control mismatch duplex. 
Samples were treated with 2% formaldehyde for lOmin 
for cross-linking. No cross-link samples were also 
prepared as negative controls. Nuclei were harvested 
using HL buffer [10 mM Tris-HCl (pH 7.5), 10 mM 
NaCl, 3mM MgCl 2 , 0.5% (v/v) NP-40]. Three million 
nuclei were resuspended in 500 ul of 1. 2 x restriction 
enzyme buffer containing 0.1% (w/v) SDS and incubated 
at 65°C for 8 min to loosen chromatin. To sequester SDS, 
Triton X-100 was added up to 1% (v/v). Samples were 
treated with 250 U of Dpnll (NEB) at 37°C overnight. 
The next day, SDS was added to 1.6% (v/v) and heated 
at 65°C for 20 min to inactivate the restriction enzyme 
followed by adding 1% (v/v) Triton X-100. For ligation, 
samples were incubated with 4000 cohesive end units of T4 
DNA ligase (NEB) at 16°C for 5h. Samples were then 
treated with 200 mM NaCl and 300 ug of proteinase 
K at 65° C overnight for reverse cross-linking. After 
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treating samples with 40ng/ul RNase A (300 jug) at 
37°C for 40min, DNA fragments were purified by 
phenol-chloroform extraction and ethanol precipitation. 
The 3C products were further purified using Wizard SV 
Gel and PCR clean-up system (Promega). Purified 3C 
samples were diluted with nuclease-free water and used 
for qPCR. iTaq SYBR Green (Biorad) was used for 
constant primer CI or C3 and PrimeTime qPCR Assay 
(IDT) was used for constant primer C2 in 3C-qPCR. 
Primers and a probe used in 3C experiments are listed in 
Supplementary Table SI. Digestion efficiency evaluated 
by qPCR amplifying across the restriction sites was 
~90%. All 3C-qPCR products were confirmed by 
sequencing. 

RESULTS 

Sense and antisense promoter RNAs are transcribed 
across the COX-2 locus 

We sequenced polyA-selected RNA from the nuclei of 
A549 non-small cell lung cancer cells to identify candidate 
genes that express long RNAs overlapping their pro- 
moters. We obtained 76 million paired-end (2 x 100) 
reads that were uniquely aligned to human genome- 19. 
We observed reads spanning the COX-2 promoter 
(Supplementary Figure SI A) and then investigated tran- 
scription at the COX-2 locus in more detail. 

We defined the origin and orientation of the transcripts 
using 5' and 3' RACE (Supplementary Figure SIB) 
and observed divergent sense and antisense domains of 
transcription at the COX-2 promoter (Figure IB, 
Supplementary Figure SIC). The COX-2 promoter is 
controlled by a TATA box and mRNA transcription ini- 
tiates almost entirely at the annotated +1 site. To quantify 
RNA levels, we used strand- specific RT-qPCR and found 
that sense transcription was prevalent near the transcrip- 
tion start site, and antisense transcription was predomin- 
ant upstream (Figure 1C). 

Both droplet digital PCR and qPCR showed that the 
sense transcripts upstream of the transcription start site 
were present in approximately two copies per cell (22). 
Transcript distribution was similar in cell nuclei and cyto- 
plasm. It is possible that the sense transcripts encode 
variants that can be translated into COX-2 protein. 
Protein produced by upstream variant transcripts, 
however, would be low relative to protein produced 
from mRNA that initiates at the +1 transcription start 
site. These data reveal an RNA network overlapping key 
regulatory regions in the COX-2 promoter. 

Coding and non-coding transcription at the COX-2 locus 
is coordinated 

The presence of an RNA network at the COX-2 promoter 
suggested a possible role for promoter RNA to regulate 
COX-2 gene expression. To investigate this possibility, we 
measured promoter RNA and mRNA expression in A549 
lung cancer cells and in the breast cancer cell lines, MDA- 
MB-231, MCF7 and T47D. These cell lines were chosen 
because their basal levels of COX-2 expression differ, 
ranging from moderate to very low. We used primer sets 



to detect RNA levels within the upstream domain of anti- 
sense transcription (primer set —535/— 387), a middle 
region where less transcription occurs (primer set —255/ 
— 193) and the downstream domain of transcription 
(primer set —81/+ 5 6). 

Expression of COX-2 mRNA and promoter RNA 
varied similarly to one another in the cell lines and 
human tissues (Figure ID and E). COX-2 expression in 
myometrial tissue is elevated just before parturition and 
provides a classic example of acute COX-2 upregulation in 
response to physiologic stress (23). In patient samples, 
increased expression of both sense and antisense RNAs 
was associated with increased expression of COX-2 
mRNA (Figure IF). 

The COX-2 promoter is responsive to endogenous 
signals including pro-inflammatory cytokines, growth 
factors and hormones. Treatment of A549 lung cancer 
cells with the COX-2 gene activators, PMA, IL1(3 and 
TNFoc increased levels of COX-2 protein (Figure 1G), 
mRNA and promoter RNA (Figure 1H). Combined, 
these results from stimuli-treated cells, cancer cells, 
patient samples and human tissues all support the conclu- 
sion that expression of COX-2 mRNA and promoter 
RNA domains vary in a coordinated and potentially 
co-regulated fashion. 

Endogenous nuclear miR-589 controls expression 
of COX-2 

It is not clear how nuclear non-coding RNAs associated 
with promoters might alter gene expression. We 
hypothesized that nuclear miRNAs could recognize 
promoter RNA and affect transcription of COX-2. To 
address this hypothesis, we first performed deep 
sequencing of nuclear RNA from A549 cells to reveal 
the identity of highly expressed miRNAs. We then per- 
formed computational searches to identify miRNAs that 
were expressed in A549 cells and possessed seed sequence 
complementarity to RNA at the COX-2 promoter. 

miR-589 was a top candidate because one strand of 
miR-589, miR-589-5p, possessed seed sequence comple- 
mentarity to two adjacent sequences within the sense 
promoter RNA (Figure 2 A and B). Both putative target 
sites are perfectly conserved in primates (Figure 2C). We 
obtained anti-miRs, short single-stranded LNA oligo- 
nucleotides, to investigate the endogenous function of 
miR-589 on COX-2 gene expression. Treatment with the 
anti-miR complementary to miR-589-5p resulted in 
reduced basal expression of COX-2 (Figure 3A). This 
result suggested that miR-589-5p may be an endogenous 
activator of COX-2. By contrast, an anti-miR complemen- 
tary to miR-452 (an miRNA with just one putative 
binding site within the COX-2 promoter, Figure 2A) had 
no effect on basal COX-2 expression. 

We then examined how introducing miR-589 into A549 
cells would affect COX-2 expression. When miR-589 was 
added as either mature miRNA with two partially com- 
plementary strands or as a miRNA mimic with two fully 
complementary strands, we observed increased levels of 
COX-2 protein (Figure 3B) and mRNA (Figure 3C). 
This increased expression of COX-2 further characterizes 
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miR-589 as an endogenous activating factor in the 
nucleus. 

Although fully complementary siRNAs cause target 
transcript levels to decrease, miR-589 has central 
mismatches relative to its target sites within the COX-2 
promoter RNA, which should prevent transcript cleavage 
by RNAi factors (24). As predicted, treatment with miR- 
589 did not decrease levels of the promoter RNA. Instead, 
promoter RNA levels increased by 1.5-3-fold (Figure 3D). 
The increase in promoter transcript levels may reflect a 
general increase in transcription at the COX-2 locus that 
is not limited to increased production of mRNA. 
Recruitment of RNA polymerase 2 (RNAP2) to the 
COX-2 promoter also increased (Figure 3E) in nuclei har- 
vested from cells treated with miR-589, consistent with 
miR-589 activating COX-2 transcription. 

To examine the action of miR-589 in a different cellular 
context, we introduced miR-589 into T47D breast cancer 
cells. Unlike A549 cells, expression of miR-589 is un- 
detectable, and COX-2 expression is substantially lower. 
We found that addition of miR-589 increased COX-2 
mRNA levels ~4-fold (Figure 3F). These data show that 



RNA-mediated activation of COX-2 expression can occur 
over a wide range of basal expression levels. 

We next sought to determine whether promoter RNA 
levels were directly linked to COX-2 expression by 
reducing expression of the promoter RNA using a 
gapmer oligonucleotide. Gapmers are antisense oligo- 
nucleotides with chemically modified flanking regions to 
increase affinity and a central DNA region to recruit 
RNase H for efficient cleavage of DNA- RNA hybrids 
(25). Gapmers are single-stranded and can selectively 
reduce the expression of either sense or antisense 
promoter RNAs. 

We identified a gapmer oligonucleotide (gapl) that 
reduced levels of the sense promoter RNA, and a 
gapmer complementary to the promoter antisense strand 
(gap2) that did not (Figure 3G). When gapl reduced levels 
of the sense promoter RNA, expression of COX-2 mRNA 
was also reduced (Figure 3H). Alterations in COX-2 ex- 
pression following treatment with anti-miR, miRNA 
mimic or gapmer are all consistent with the sense 
promoter RNA acting as a docking site for miR-589. 
Both promoter RNA and miR-589 are required for 
RNA-mediated activation of COX-2 expression. 
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Transcriptional activation of the COX-2 gene by miR-589 
involves the RNAi factors, AG02 and GW182 

Since small RNAs require RNAi proteins to function, we 
reasoned that the RNAi factor argonaute-2 (AG02) 



protein might be involved in RNA-mediated activation 
of COX-2. It is well known that AG02 promotes 
binding of small RNAs to RNA transcripts during cyto- 
plasmic RNAi. AG02 is also found in cell nuclei of yeast, 
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plants and animals (26). RNAi has been reported to 
function in mammalian nuclei (27) where AG02 can 
facilitate RNA-mediated modulation of transcription 
(20,21) or splicing (28). 

Using stringent cellular fractionation protocols 
followed by western blot analysis, we found that both 
AG02 and a second RNAi factor, GW182 (29), are 
present in the nuclei of A549 cells (Figure 4A). We per- 
formed RIP from nuclear extracts using an anti-AG02 
antibody to determine whether AG02 associates with 
promoter RNA at the COX-2 promoter. RIP revealed 
the recruitment of AG02 to the COX-2 promoter RNA 
was dependent on the addition of miR-589 (Figure 4B). 

To further investigate a potential role for AG02 in 
transcriptional activation by small RNA, we used anti- 
AG02 siRNAs to deplete AG02 levels inside cells and 
examined how reduction of AG02 expression would 
affect gene activation. We observed that RNA-mediated 
gene activation by miR-589 was blocked when AG02 
levels were reduced (Figure 4C). Similar experiments 
using an RNA with seed sequence mismatches did not 



block gene activation or cause AG02 recruitment. These 
two complementary experimental approaches, RIP and 
AG02 depletion, support the conclusion that AG02 is 
required for the sequence-selective recognition of miR- 
589 to the promoter RNA and support a role for 
nuclear RNAi factors in endogenous transcriptional 
regulation. 

GW182 (TNRC6A) is an important factor in the mech- 
anism of translational silencing by miRNAs (29). GW182 
interacts directly with AG02 through proteimprotein 
interactions (30-32). Co-immunoprecipitations with anti- 
AG02 or anti-GW182 antibody revealed that AG02 and 
GW182 form a complex in stringently purified cell nuclei 
(Figure 4D). To determine whether activation of COX-2 
expression requires GW182, we used an siRNA pool 
capable of reducing expression of all three GW182 
paralogs (Supplementary Table S2). Reducing GW182 ex- 
pression reversed miR-589-mediated activation of COX-2 
expression relative to A549 cells that had been treated 
with a control RNA that left GW182 levels unchanged 
(Figure 4E). 
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These results indicate that the RNAi factors AG02 and 
GW182 associate in the nucleus (Figure 4D and F) and are 
required for transcriptional activation of the COX-2 gene 
by miR-589 (Figure 4C and E). Taken together, the results 
suggest that AG02 forms a complex with GW182 and 
miR-589 that is guided to the promoter RNA. Binding 
specificity is achieved through Watson-Crick base- 
pairing of miR-589 to the promoter RNA (Figure 4F). 

Discovery of more potent activating small RNAs 

Activation of COX-2 gene expression by miR-589 sug- 
gested that it would be possible to discover related small 
RNAs that were more potent modulators of gene expres- 
sion. Selective and robust RNA activators would be trans- 
formative as: (i) tools to investigate the mechanism of 
RNA-mediated gene activation (higher levels of activation 
increase signal to noise and facilitate the diverse experi- 
ments); (ii) facilitators for gain of function studies in the 
same manner that RNAi knock-down enables loss of 
function experiments; and (iii) new options in the devel- 
opment of therapeutics where upregulation of a specific 
gene is desired. Duplex RNAs that are fully complemen- 
tary to gene promoters have been reported to activate 
several genes, suggesting that modulation of transcription 
by small RNAs may be a general regulatory mechanism 
(5,20,21,33-38). 

We designed small RNAs complementary to sequences 
adjacent to or overlapping the target sites of miR-589 
within the COX-2 promoter RNA (Figure 5A). In 
contrast to miR-589, these duplex RNAs had a single 
target site and were fully complementary. We found that 
several of these small RNAs activated COX-2 expression 
(Figure 5B). 

One of these RNAs, RNA 12, achieved 25-fold activa- 
tion, significantly more than the 8-fold enhancement 
achieved by miR-589. Activation was potent with half 
maximal activation achieved at < 10 nM (Figure 5C). The 
seed sequence necessary for efficient recognition by small 
RNAs during RNAi occupies bases 2-8 within the guide 
strand of a duplex RNA (39). The introduction of one or 
more mismatched bases into the seed sequence of RNA 12 
dramatically reduced activity (Figure 5D), consistent with 
the promoter RNA being the molecular target. Like miR- 
589, RNA 12 activated COX-2 expression in T47D cells 
(Figure 5E). Combined treatment of RNA 12 and miR- 
589 did not result in additive activation of COX-2 expres- 
sion (Supplementary Figure S2A). Activation was transi- 
ent, decreasing after 7 days (Figure 5F). Although 
transient, the activation was also persistent, as it 
remained detectable after 9 days. 

Although RNA 12 was the strongest activator of COX-2 
expression, it was not the only RNA activator discovered 
in our screen. Several similar RNAs also activate expres- 
sion by as much as 5-10-fold (Figure 5B, RNA10, 
RNAI 3, RNA38). One of these RNAs, RNA38, possesses 
a seed sequence that overlaps with miR-589. Other RNAs 
were less active, consistent with previous observations that 
closely related nucleic acid sequences can vary in activity 
(40). This variation may be due to differences in the ac- 
cessibility of the target RNA, association or usage by 



RNAi factors or the ability to form productive inter- 
actions at the promoter. We have previously observed 
that, when introduced into cells first, inactive duplex 
RNAs that are fully complementary to a promoter tran- 
script blocked gene upregulation of progesterone receptor 
by an activating RNA with an overlapping target site (34). 
This result shows that the exact location for recognition 
matters during activation. 

The large dynamic range in gene activation by RNA 12 
made it an ideal small RNA for evaluating mechanism. 
RNA 12 treatment increased COX-2 pre-mRNA levels, a 
result consistent with increased gene transcription 
(Supplementary Figure S2B). ChIP from purified nu- 
clei revealed that activation by RNA 12 induced recruit- 
ment of RNAP2 (Figure 5G), CREB1 and NFkB 
(p65, p50) (Supplementary Figure S3) to the COX-2 
promoter. Recruitment of NFkB and CREB1 on activa- 
tion of the COX-2 gene is also consistent with induction of 
transcription by RNA 12 (16). Taken together, these 
observations provide support for the existence of an 
RNA-mediated mechanism that can activate gene 
transcription. 

When RNA strands are introduced into cells, they can 
provoke an interferon response and that might induce 
COX-2 expression. To address this possibility, we 
assayed the activity of an RNA 12 variant (RNA 12m 
m = modified). Like RNA 12, RNA 12m is fully comple- 
mentary to its target but differs by introduction of 
chemical modifications known to minimize the interferon 
response (41). RNA 12m activated COX-2 expression by 
~ 10-fold (Supplementary Figure S2C), demonstrating 
that RNA designs compatible with in vivo development 
are effective activators of gene expression. Levels of inter- 
feron responsive genes after treatment with RNA 12 or 
RNA 12m did not significantly increase (Supplementary 
Figure S4A and B). These data, together with stringent 
discrimination against activation by duplexes with seed 
sequence mismatches, strongly support the hypothesis 
that small RNAs targeted to the COX-2 transcription 
start site are ' on- target' activators of COX-2 gene 
expression. 

Small RNAs mediate a distinct pathway for COX-2 gene 
activation 

COX-2 expression can be activated by the known pro-in- 
flammatory factors IL1(3, TNFoc and PMA, or repressed 
by the anti-inflammatory factor, Dexamethasone (Dex). 
To understand how the mechanism of RNA-mediated ac- 
tivation compares with these pathways, we treated A549 
cells with these compounds alone or in combination with 
RNA12. 

IL1(3 increased COX-2 expression to a maximum level 
of 15-20-fold (Supplementary Figure S2D, Figure 6A), 
increased COX-2 pre-mRNA levels (Supplementary 
Figure S2B), and increased recruitment of RNAP2 at 
the COX-2 promoter (Figure 5G). TNFoc and PMA 
increased COX-2 levels by 10- and 8-fold, respectively 
(Figure 6A). Compared with treatment with IL1(3, 
TNFoc or PMA alone, when cells were treated cells with 
RNA 12 in combination with IL1(3, TNFoc or PMA, we 
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observed increased COX-2 expression to higher levels, 
106-, 93- and 53-fold, respectively (Figure 6A). Co-treat- 
ment with RNA 12, therefore, yields expression levels that 
are much greater than can be achieved by the maximal 
doses of known activators of COX-2. 

We also tested RNA 12 in combination with Dex an 
anti-inflammatory agent known to repress COX-2 expres- 
sion. Addition of Dex reversed the effects of known acti- 
vators PMA, ILip and TNFoc (Figure 6B). RNA 12, by 
contrast, overrides repression by Dex. 

Super-activation of COX-2 expression beyond what can 
be induced by known activators, combined with the unique 
ability to override repression by Dex, indicates that the 
mechanism of promoter RNA-mediated activation is 
distinct from the known mechanisms of other inflammatory 
pathway regulators. Small RNAs provide a new sequence- 
specific pathway for activating COX-2 expression that 
requires complementarity to the COX-2 promoter RNA. 



Sense promoter RNA is the molecular target for COX-2 
gene activation 

Only one strand of miR-589 is complementary to the 
promoter RNA, making the sense strand the only poten- 
tial partner for gene activation. RNA 12, however, is com- 
plementary to both sense and antisense promoter RNAs 
making it possible that either stand might be its target. 

To determine whether RNA 12, like miR-589, targets 
the sense strand, we mutated a single base within the 
seed sequence on either strand of RNA 12. We annealed 
the strands containing a mutation with native complemen- 
tary strands lacking a mutation (Figure 7A). When the 
strand complementary to the sense promoter RNA was 
mutated at position 3 in the seed sequence, gene activation 
was almost abolished. Mutation of the strand complemen- 
tary to the antisense strand had little impact. These data 
suggest that the sense strand of the promoter RNA is the 
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molecular target for RNA-mediated activation of COX-2 
expression. 

We also used antisense oligonucleotide gapl to reduce 
expression of the sense RNA target and then treated cells 
with RNA 12. Antisense gapmer oligonucleotides reduce 
promoter sense transcripts by inducing degradation of 
RNA-DNA targets through recruitment of RNase H. 
Because gapmers are single-stranded, they are good tools 
for distinguishing differences between the sense and anti- 
sense transcripts. 

When RNA 12 is introduced following transfection of 
A549 cells with gapl, we observed decreased activation 
of COX-2 expression from 7- to 2.5-fold (Figure 7B). 
This result is consistent with promoter RNA acting as a 
scaffold that is necessary for small RNA-mediated mo- 
lecular recognition and subsequent activation of COX-2 
expression. By reducing sense transcript levels, the gapmer 
removes many of the target sense transcripts for subse- 
quent recognition by AG02/small RNA complexes and 
would therefore block gene activation. 

Taken together, three lines of evidence support the con- 
clusion that COX-2 gene activation by either RNA 12 or 
miR-589 involves recognition of the sense promoter RNA, 
(i) miR-589 is complementary to the sense promoter tran- 
script, not the antisense transcript; (ii) a single mutation of 
the seed sequence within the antisense strand of RNA 12 
blocks gene activation; and (iii) antisense oligonucleotide- 
mediated reduction of sense strand expression reverses 
gene activation. 

Promoter RNA acts as a scaffold for small RNA/RNAi 
factor complexes 



To further explore the potential for mechanistic similarity 
between miR-589 and RNA 12, we examined the role of 
RNAi factors. Similar to miR-589, activation by RNA 12 
also involved RNAi factors. RIP using nuclear extracts 
revealed that RNA 12 recruited AG02 to COX-2 
promoter RNA (Figure 8 A). When AG02 levels 
are reduced using an anti-AG02 siRNA, promoter 



RNA-mediated gene activation by RNA 12 is blocked 
(Figure 8B). When GW182 levels are reduced using a 
siRNA pool to silence all three TNRC6 paralogs, 
promoter RNA-mediated gene activation by RNA 12 
is abolished (Figure 8C, Supplementary Figure S5). 
This requirement for AG02 and GW182 parallels 
the mechanism of miR-589 action described earlier in 
the text. 

Cleavage of the promoter RNA is not necessary for gene 
activation 

miR-589 and RNA12 share a requirement for RNAi 
factors. RNA 12, however, is fully complementary and 
would be expected to induce AG02-mediated cleavage 
of the promoter RNA (42). Consistent with this expect- 
ation, we observed that transfection of RNA 12 into A549 
cells reduces promoter RNA levels (Figure 9A). miR-589, 
by contrast, contains mismatches that disrupt the poten- 
tial for cleavage by AG02 (24) and does not reduce 
promoter RNA levels (Figure 3D). 

To further investigate whether target cleavage was 
involved or required in the gene activation mechanism, 
we made RNA 12 more 'miR-like' by introducing central 
mismatches to produce RNA12^c (nc = no cleavage) 
(Figure 9B). The two synthetic RNAs (RNA 12 and 
RNA12^c) were transfected into cells. Nuclear or cyto- 
plasmic extracts were prepared and used to test cleavage 
of a radiolabeled substrate RNA containing the common 
target site for RNA 12 and RNAllnc. 

RNA 12 induced cleavage of the RNA substrate, 
whereas KNAllnc did not. Consistent with its failure to 
cleave a synthetic substrate and in contrast to RNA 12 
(Figure 9 A), KNAllnc did not reduce promoter RNA 
levels (Figure 9C). Instead, similar to the effect on 
promoter RNA on addition of miR-589 (Figure 3D), 
addition of KNAllnc increased promoter RNA levels 
5-fold. 

RNA 12 and RNA12^c have opposing effects on the 
level of promoter RNA, but their effects on COX-2 
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Figure 8. Gene activation by fully complementary RNA involves recruitment of RNAi factors. (A) RIP showing recruitment of AG02 protein to 
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expression and their mechanism of action appear to be the 
same. When added to A549 or T47D cells, RNAUnc 
activated COX-2 expression as efficiently as RNA 12 
(Figures 5E and 9D). Like RNA 12, RNAllnc super- 
activated COX-2 expression when combined with IL1(3, 
TNF or PMA (Supplementary Figure S2E). RIP using 
nuclear extracts demonstrated recruitment of GW182 to 
the promoter RNA (Figure 9E). 

Our results showing that KNAllnc can induce activa- 
tion of COX-2 expression reinforce the data from miR-589 
showing that activation can occur in the absence of target 
cleavage. These results further support the conclusion 
that RNA at the COX-2 promoter acts as a scaffold for 
binding the small RNA/RNAi factor complexes 
that trigger gene activation. After induction of this gene 
activation cascade, the potential for cleavage of the 
promoter RNA and its absolute levels are not critical to 
mechanism. 



Gene activation involves modification of histones via 
WDR5 

To test the effects of promoter RNA-mediated gene acti- 
vation on histone modification, we performed ChIP using 
primer sets covering regions from the +1 transcription 
start site to —1540 nucleotides upstream (Supplementary 
Figure S6A). H3K4 trimethylation and H4 acetylation 
are associated with gene activation. ChIP analysis using 
purified nuclear extract revealed both histone marks were 
increased following treatment with RNA 12 or KNAllnc 
(Figure 10A, Supplementary Figure S6B and C). The oc- 
currence of repressive mark H3K27 trimethylation was 
unchanged (Figure 10A, Supplementary Figure S6D). 

WDR5 is a WD40 repeat-containing protein that can act 
as a protein scaffold to stimulate histone methyltransferase 
activity (43). WDR5 has been reported to associate with 
HOTTIP, an IncRNA that controls activation of several 
HOXA genes (44). To determine whether WDR5 was 
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involved in RNA-mediated control of COX-2 gene expres- 
sion, we used an siRNA against WDR5 mRNA to decrease 
its expression (Supplementary Figure S6E). 

We found that reduced levels of WDR5 blocked RNA- 
mediated gene activation of both COX-2 protein and 
mRNA expression by RNA 12 and RNA12ftc (Figure 10B 
and C). Reduced levels of WDR5 lowered the induction of 
H3K4me3 modification of histones by RNA12 and KNAUnc 
in purified nuclear extracts (Figure 10D). Also consistent with 
RNA-mediated action of WDR5, ChIP from nuclear extracts 
indicated that addition of RNA 12 increased levels of WDR5 
recruitment at the COX-2 promoter (Figure 10E). These 
results suggest that RNA-mediated gene activation by small 
RNA/RNAi factor complexes induces chromatin modifica- 
tion and requires WDR5. 

Promoter RNA links activation of COX-2 and PLA2G4A 
gene transcription 

COX-2 is one of several key enzymes in the eicosanoid 
pathway that act in concert to regulate inflammation. 



PLA2G4A catalyzes the hydrolysis of membrane 
phospholipids to release arachidonic acid, the physiologic 
substrate of COX-2 (Figure 11 A). The gene 
encoding PLA2G4A protein is adjacent to the COX-2 
gene on chromosome 1 and is M49kb distant, 
indicating that the close functional partnership be- 
tween the two genes may be accompanied by spatial 
proximity. 

The physiologic connection and proximal relationship 
between COX-2 and PLA2G4A led us to test the hypoth- 
esis that activating RNAs might also control PLA2G4A 
expression. Addition of the anti-miR targeted to miR-589 
decreased basal PLA2G4A expression (Figure 1 IB). 
Conversely, addition of miR-589 led to an increase in 
basal expression levels (Figure 11C). RNA 12 and 
RNA12ftc also activated PLA2G4A expression, once 
again showing their functional equivalence (Figure 11D). 
As observed for activation of COX-2 expression, insertion 
of seed sequence mismatches abolished activation of 
PLA2G4A expression. ChIP analysis from nuclear 
extracts revealed increased recruitment of RNAP2 at the 
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PLA2G4A promoter, a result consistent with transcrip- 
tional activation by small RNA (Figure HE). 

Two conclusions can be drawn from these results. The 
first is that small RNAs that sequence-selectively activate 
COX-2 expression also activate expression of PLA2G4A, 
the gene responsible for the enzyme that yields the sub- 
strate for the COX-2 enzyme. Activation of PLA2G4A 
mirrored activation of COX-2. The second conclusion is 
that a common mechanism of RNA-mediated regulation 
exists for the two genes. 

Similarities in the mechanism of activation for 
PLA2G4A and COX-2 include the following: (i) addition 
of antisense oligonucleotide gapl to reduce levels of the 
sense promoter RNA reversed PLA2G4A activation by 
RNA 12 (Figure 1 IF); (ii) siRNA-mediated reduction of 
AG02, GW182 or WDR5 reversed activation of 
PLA2G4A expression (Figures 11GHI, Supplementary 
Figure S7A); (iii) addition of small RNA increased levels 
of H3K4me3 modifications throughout the PLA2G4A 
promoter, consistent with involvement of WDR5 in 
RNA-mediated activation (Supplementary Figure S7B 
and C); and (v) activation increased H4Ac marks across 
the PLA2G4A promoter (Supplementary Figure S7B and 
D). These mechanistic similarities suggest a link between 
the two genes that can be exploited by RNA. 

It was possible that increased levels of COX-2 protein 
or enzymatic activity might explain increased expression 



of PLA2G4A. To test this hypothesis, we treated A549 
cells with activators IL1(3, TNFoc or PMA to elicit 
maximal activation of COX-2 protein expression and 
observed a less than 2-fold increase in expression of 
PLA2G4A (Supplementary Figure S7E). Conversely, 
blocking COX-2 enzymatic activity by addition of inhibi- 
tor did not alter RNA-mediated activation of COX-2 
(Supplementary Figure S8A) or PLA2G4A 
(Supplementary Figure S8B). 

Physical linkage between COX-2 and PLA2G4A 
promoters 

One explanation for increased expression of both COX-2 
and PLA2G4A would be a physical linkage between the 
two genes. Gene looping can occur within a gene, between 
a gene and an enhancer site or between genes (45) to fa- 
cilitate the action of noncoding RNAs (36). The COX-2 
and PLA2G4A gene loci map to chromosome lq25 in a 
head to head configuration separated by 149 kb. There is 
promoter RNA overlapping the COX-2 promoter, but no 
RNA was detected overlapping the PLA2G4A promoter 
in the RNA-Seq analysis. Gene looping would explain 
why interactions at the COX-2 promoter RNA might 
affect PLA2G4A. 

To determine whether COX-2 and PLA2G4A loci 
associate, we performed chromosome conformation 
capture (3C) analysis (46). In 3C, chromosomal DNA is 
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Figure 11. Promoter RNAs link COX-2 and PLA2G4A transcription. (A) Distal and functional relationship of COX-2 and PLA2G4A. Scheme 
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cross-linked and treated with a restriction enzyme. 
Sequences that are held in close proximity by chr- 
omosomal structure, but which might be far distant 
from one another on the linear chromosome, can be 
ligated and subsequently detected by PCR and DNA 
sequencing. 

We performed 3C analysis with one primer (CI) at the 
COX-2 promoter and other primers targeting the 
PLA2G4A promoter (primers T4 or T5) (Figure 12A, 
Supplementary Table SI). All experiments used purified 
A549 nuclei. Using SYBR-based PCR, a method that 
allows for the qualitative detection of cross-linked 
products, we observed amplified products in samples 



that had been cross-linked, but not in samples that were 
not cross-linked (Figure 12B). 

The S'-nuclease assays including PrimeTime (IDT) and 
TaqMan (ABI) qPCR are more suitable than the SYBR- 
based qPCR for 3C because a target-specific fluorescent 
probe for detection supplements gene-specific primers 
used to initiate product amplification. Use of a fluorescent 
probe for detection reduces the likelihood of detecting 
non-specific products and makes quantitative determin- 
ations more accurate. Using primers (C2 and Tl-8) tar- 
geting the regions between COX-2 and PLA2G4A genes 
and a probe for qPCR (Figure 12A), we determined 
whether associations were being selectively enriched 
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near the transcription start sites for COX-2 and 
PLA2G4A. 

As observed when using SYBR chemistry, cross-linked 
products could be observed in the presence, but not in the 
absence, of cross-linking (Figure 12C). Quantitative meas- 
urement of amplification revealed that RNA12 selectively 
increased association between the COX-2 and PLA2G4A 
promoters (Figure 12D). These results are consistent with 
the suggestion that RNA-mediated activation of COX-2 
and PLA2G4A by RNA12 is associated with closer prox- 
imity between the two promoters. 

Sequencing the isolated ligation products confirmed 
junctions between the COX-2 and PLA2G4A promoters 
(Figure 12E), indicating that the two regions were in close 
proximity in chromosomal DNA. The close physical prox- 
imity of the COX-2 and PLA2G4A promoters provides a 
mechanism for communicating activation signals over 
long genomic distances and allows the functional partner- 
ship of the two enzymes during eicosanoid biosynthesis to 
be accompanied by linked expression. 

PLA2G4A is just one of several PLA2 family genes. We 
examined whether the COX-2 promoter might also be 
physically linked to these related genes. 3C analysis with 
multiple primer sets (CI or C3; T9-12) showed no 
evidence for connections between the COX-2 promoter 
and genes encoding extracellular PLA2 that cluster on 
the opposite arm of chromosome 1, ~166 million bases 
distant (Figures 12FG). 

DISCUSSION 

RNAi factors interact with RNA in cell nuclei to activate 
gene expression 

miRNAs and RNAi factors are present in mammalian 
cell nuclei. These factors can form interactions with each 
other and with RNAs that overlap gene promoters 
(Figure 13 A). Our results are best explained by a mechan- 
ism involving the RNA at the COX-2 promoter acting as a 
scaffold for small RNA-mediated recruitment of AG02 
and GW182 near the COX-2 promoter in cis 
(Figure 13B). Similar involvement of RNAi factors and 
nascent RNA scaffolds during transcriptional modulation 
and heterochromatin assembly in yeast has been noted 
(47,48). Recognition of the promoter RNA by the RNA- 
AG02-GW182 complex may trigger recruitment of other 
factors like WDR5 and activation of COX-2 and 
PLA2G4A gene expression. The 3C analysis reveals inter- 
action between the COX-2 and PLA2G4A promoters, 
demonstrating a gene architecture that permits long 
distance signaling (Figure 13C). 



Activation of COX-2 expression requires on-target 
recognition at the COX-2 loci 

Whenever nucleic acids are introduced into cells there is a 
potential to produce off-target effects. These effects can 
confound experiments by causing misleading changes in 
target gene expression. Guarding against off-target 
effects is essential for any experiment with duplex RNA. 

The use of multiple mismatched or scrambled control 
duplexes is a convenient and effective method to establish 
whether gene activation requires complementarity. We 
find that as little as one mismatch within the seed 
sequence largely abolishes gene activation (Figure 7A), 
consistent with an on-target effect that functions 
through the RNAi pathway. Scrambled sequences or se- 
quences with multiple mismatches spread throughout the 
RNA were also inactive. 

When RNA is introduced into cells, it can induce an 
interferon response that might be especially problematic 
for a highly inducible gene like COX-2. We observe little 
or no increase in the expression of a panel of interferon 
responsive genes (Supplementary Figure S4B). 

COX-2 gene activation is also achieved in the presence 
of RNA duplexes that contain modified nucleosides that 
are designed to reduce off-target effects (Supplementary 
Figure S2C). Addition of siRNAs designed to reduce 
AG02 or GW182 expression, or a gapmer oligonucleotide 
designed to reduce levels of the target promoter transcript, 
reduce COX-2 expression rather than increase it as would 
be expected from an off-target effect that occurs after 
addition of exogenous nucleic acid. 

RNA12/RNA12rcc and miR-589 have adjacent target 
sites near the COX-2 transcription start site (Figure 5A), 
consistent with their similar ability to activate COX-2 
gene expression. These targets sites, however, have 
different seed sequence matches. The fact that recogni- 
tion by duplex RNAs with different seed sequence 
complementarities causes the same effect on expression 
of both PLA2G4A and COX-2 further supports the con- 
clusion that the effect is on-target. RIP using purified 
nuclei establishes that both RNA12/RNA12ftc and miR- 
589 associate with the target COX-2 promoter RNA, sup- 
porting an on-target physical interaction. RNA 12 and 
RNA12/7C have different sequences and different potential 
to induce cleavage of complementary target sites. The 
finding that their effects on COX-2 expression are virtually 
identical further supports the conclusion that the observed 
activation is on-target. 

The details of RNA-mediated gene activation differ 
from activation by known activators like ILip, TNFoc 
and PMA (Figure 6A). These data rule out an off-target 
effect through known pathways for COX-2 activation. 
There is no known direct link between COX-2 and 



Figure 12. Continued 

(7-test) relative to MM. (E) DNA sequencing of amplified 3C PCR products (C1+T4, C1+T5, C2 + T5) demonstrates a junction between the 
PLA2G4A and COX-2 promoter sequences. (F) Diagram of COX-2/PLA2G2E/PLA2G2A gene loci and design of 3C experiments. The 3C PCRs 
were performed using combinations of constant primers (CI or C3) and test primers [T9, T10 (PLA2G2E promoter); Til, T12 (PLA2G2A 
promoter)]. Mismatch control or RNA12 were transfected at 25 nM. (G) Analysis of 3C PCR products on 2.5% agarose gels. Expected amplicon 
sizes for each primer set were 168-bp (C1 + T9), 131-bp (C1+T10), 142-bp (C3 + T9), 105-bp (C3 + T10), 167-bp (Cl + Tll), 165-bp (C1 + T12), 
141-bp (C3 + T11) and 139-bp (C3 + T12), but no specific products were detected after 40 cycles of PCR amplification. 
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Figure 13. Model of promoter RNA-mediated activation of COX-2 and PLA2G4A genes. (A) Ground state: promoter RNAs are expressed at the 
COX-2 promoter. COX-2 and PLA2G4A are expressed at low levels. (B) Intermediate state: The sense promoter RNA is a platform for recognition 
by a small complementary RNA (endogenous or synthetic) and recruitment of AG02 and GW182. (C) Gene activation state: Binding of small RNA/ 
AG02/GW182 complex leads to association of transcription factors and WDR5, and an increase in activating histone marks. Specific gene looping is 
induced between the COX-2 and PLA2G4A promoters, enabling activations of both genes. 
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PLA2G4A expression. Any off-target cause for linked 
upregulation of COX-2 and PLA2G4A would be unprece- 
dented. Taken together, our data and the robustness of 
RNA-mediated gene activation of COX-2 and PLA2G4A 
make a strong case for on-target gene activation. 

Scaffold recognition, not strand cleavage, triggers 
activation 

Although miR-589, RNA12 and RNA12^c are comple- 
mentary to adjacent target sites, the three RNAs have 
fundamental differences in their structures and potential 
for affecting gene expression. They bind different se- 
quences and differ in their complementarity relative to 
the promoter RNA. By evaluating their differences and 
similarities, it is possible to draw several conclusions 
about mechanism. 

miR-589 and RNA12^c contain central mismatched 
bases that prevent cleavage of target transcripts by 
AG02. RNA 12, by contrast, is fully complementary and 
enables AG02 to cleave the target RNA (Figure 9B). Our 
experimental observations are consistent with predicted 
potentials for strand cleavage. RNA 12 reduces levels of 
COX-2 promoter RNA (Figure 9A) and induces 
cleavage of target RNA (Figure 9B), as would be 
expected for a fully complementary RNA/AG02 
complex. miR-589 and RNA12/?c, which both contain 
central mismatches, do not reduce levels of the promoter 
RNA (Figure 3D, Figure 9C) or induce cleavage 
(Figure 9B). 

The potential for cleaving the target transcript is a fun- 
damental difference between RNA 12 and miR-589/ 
RNA12fzc. Despite this fundamental difference, however, 
their activities are strikingly similar. RNA12/RNA12/?c 
and miR-589 (i) activate COX-2 and PLA2G4A expres- 
sion; (ii) require seed sequence complementarity for acti- 
vation; (iii) recruit AG02 to the promoter RNA, (iv) 
induce histone modifications; (v) have their activity 
reduced by an antisense gapmer oligonucleotide that 
targets the promoter transcript; and (vi) require AG02 
and GW182 expression for activity. 

The functional similarity of miR-589, RNA 12 and 
RNA12ftc, despite their different potential to induce 
cleavage of target, demonstrates a central feature of mech- 
anism. Cleavage of the target RNA is not required for 
gene activation. Instead, the initial recognition step of 
the COX-2 promoter RNA by the small RNA and 
associated protein factors is the critical switch that 
triggers activation. The subsequent fate of the promoter 
RNA after it binds the RNA/AG02 complex, cleaved or 
uncleaved, has little effect on activation. Our results 
suggest that the most likely role for the promoter RNA 
is to act as a scaffold. 

Antisense oligonucleotide gapl reduces levels of the 
promoter transcript by RNase H-mediated cleavage. 
Reduction of the promoter transcript by gapl does not 
increase COX-2 expression (Figure 3H), further support- 
ing the conclusion that cleavage of the transcript is not 
necessary for activation. Indeed, we observe gapl leads to 
a decrease of COX-2 expression, consistent with the 



promoter transcript being needed as a scaffold during 
activation. 

Both RNA 12 and antisense oligonucleotide gapl reduce 
promoter RNA levels (Figure 9 A, Figure 3G), but only 
RNA 12 causes activation. This may seem contradictory, 
but the contradiction is resolved when the different mech- 
anisms of the two nucleic acids are considered, gapl is an 
antisense oligonucleotide, not a duplex RNA. Because it is 
a single-stranded RNA-DNA hybrid, it does not recruit 
AG02. By contrast, when RNA 12 recognizes the 
promoter transcript, it recruits AG02 and other factors 
such as GW182 (Figure 8A). 

This AG02/RNA complex recognition step, which 
gapl does not share but which is shared by RNA12/-/C 
and miR-589, is the likely trigger for gene activation. 
When gapl was added before RNA 12, the addition of 
gapl caused loss of the transcript before RNA 12 can 
recruit AG02. By removing the scaffold and the potential 
to recruit AG02, activation is prevented. Activation 
requires an intact transcript for recognition and a small 
RNA capable of recruiting AG02 and associated protein 
factors to the target. 

Promoter RNAs are targets for ribonucleoprotein 
transcription factors 

RNAs that overlap gene promoters could influence gene 
expression through several different mechanisms. These 
include the following: (i) transcriptional interference 
caused when synthesis of the IncRNA interferes with tran- 
scription of mRNA; (ii) direct hybridization between an 
antisense transcript and an mRNA; (iii) miRNAs or other 
functional RNAs encoded within IncRNAs; and (iv) 
global recruitment of histone or DNA modifications. All 
of these mechanisms may be physiologically relevant in 
mammalian cells, but our results with COX-2 / PLA2G4A 
activation through promoter RNA recognition demon- 
strate a novel mechanism — small RNA-dependent recruit- 
ment of factors at a promoter in cis to affect the protein 
machinery and the transcription of a specific gene. 

Efficient RNA recognition in the nucleus might have 
advantages as a mechanism for control of splicing or tran- 
scription. Although proteins are powerful agents for 
controlling expression, they lack a general capacity to be 
highly selective for just one gene or to rapidly evolve 
specificities for recognizing new genes. Because 
IncRNAs, miRNAs, AG02 and GW182 exist in the 
nucleus, they are the basis for an alternate regulatory 
pathway that exploits the versatility of Watson-Crick 
base-pairing. When AG02 binds a small RNA, it forms 
a ribonucleoprotein complex and facilitates recognition of 
complementary targets within nuclear RNAs. Sequence- 
specific recognition using Watson-Crick base pairing can 
evolve new specificities relatively rapidly using conserved 
core RNA/protein machinery. 

When IncRNAs are synthesized as nascent transcripts, 
they are close to gene promoters and poised to influence 
transcription in cis. Even if present at just a few copies per 
cell, the effective concentration of nascent transcripts 
relative to promoter DNA is high. Like a transcription 
factor binding to a promoter, or a ligand changing the 
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conformation of a nuclear hormone receptor, delivery of 
the RNA/protein complex to the promoter RNA scaffold 
triggers gene activation. We have recently observed that 
the COX-2 promoter transcript is present at ~2 copies per 
A549 lung cancer cell (22), consistent with action in cis but 
not consistent with the expectation that the transcript 
might locate its target in trans. 

Implications for regulation of COX-2 

It seems logical that PLA2G4A and COX-2 expression 
should be linked, as COX-2 enzyme uses the product of 
PLA2G4A enzyme. However, despite of the central role of 
eicosanoid signaling in inflammation and other critical 
physiologic processes, mechanisms for coordinating their 
expression have not been described. We find that COX-2 
and PLA2G4A are both induced by an RNA-mediated 
mechanism that can be triggered by endogenous miR- 
589. RNA provides a mechanism that allows coordinated 
expression of two key enzymes within the eicosanoid 
pathway. 

RNA-mediated induction coordinates COX-2 and 
PLA2G4A protein levels selectively, unlike the more 
global gene induction triggered by IL1(3 and other 
known inducers. For example, COX-1 and COX-2 differ 
dramatically in their regulation. COX-2 expression is 
usually induced while COX-1 is constitutively expressed. 
Gene looping and recognition of the COX-2 promoter 
RNA by miR-589/RNAi factor complex provide a 
means for coordinated expression of COX-2, but not 
COX-1, to PLA2G4A synthesis. 

Addressing an alternate hypothesis to explain COX-2/ 
PLA2G4A activation 

Any experiment that uses nucleic acids inside cells runs the 
risk of inducing off-target effects. As noted earlier in the 
text, we have performed several experiments to minimize 
this possibility. One alternate explanation, however, 
would be if miR-589 and the other five activating RNAs 
were modulating expression of one or more transcription 
factors, which in turn increased both COX-2 and 
PLA2G4A expression. 

We believe that this is unlikely for the following 
reasons: (i) As miR-589, RNA 12, RNAUnc, RNA 10, 
RNA 13 and RNA38 all activate COX-2 expression, they 
would all need to target a transcription factor that 
controls both COX-2 and PLA2G4A; (ii) None of the 
RNA activators have seed complementary to known tran- 
scriptional repressors of COX-2 or PLA2G4A; (iii) As ac- 
tivation by RNA in combination with well-established 
physiologic COX-2 activators is additive, the putative 
transcription factor would need to control a new 
pathway for COX-2 activation; (iv) An oligonucleotide 
gapmer that reduces expression of the target promoter 
transcript also reduces gene activation by miR-589 and 
RNA12. If activation was being modulated through 
downregulation of a transcription factor, targeting the 
promoter transcript would be expected to have no effect; 
and (v) We observe physical interactions between the 
promoter transcript and AG02. We also observe inter- 
actions between the promoters for PLA2G4A and 



COX-2. These contacts describe direct physical inter- 
actions between the major factors involved in our 
proposed mechanism of action, reducing the need to 
invoke an indirect mechanism that involves unidentified 
transcription factors. 

It could also be argued that the small RNAs are 
titrating a repressor that normally binds to the COX-2 
and PLA2G4A promoters. For such binding to affect 
gene expression in our experiments, it would need to be 
acutely sensitive to just a single change to one nucleoside 
of the duplex within the region normally associated with 
seed sequence activity for RNAi because one substitution 
abolishes activation (Figure 7A) or to the introduction of 
other mismatched bases (Figure 5D, Figure 11D). The 
repressor titration would need to be insensitive to intro- 
duction of two mismatched bases (activation is still 
achieved by RNA12nc) or more extensive disruption of 
base-pairing (activation by incompletely matched miR- 
589). The repressor would also need to be titrated by 
two families of RNAs (miR-589 and RNA12/RNA12rcc/ 
RNA 1 0/RNA 1 3 /RNA38) that share little overlap in their 
sequences. The repressor would need to control both 
PLA2G4A and COX-2, even though the promoters share 
no recognized repressor binding sites. Repressors 
normally bind B-form DNA, but for titration to occur 
in this case they would need to recognize a short 
RNA that can be chemically modified (Supplementary 
Figure 2C). The duplex RNA would also need to survive 
inside mammalian cells without being degraded or loaded 
into the RISC complex as single strand. Finally, a mech- 
anism of inactivation that involves titration of a repressor 
would be inconsistent with dependence on AG02 or 
GW182 or the observed reversal of activation on 
addition of a gapmer targeting the promoter transcript. 

Finally, we note that activation of COX-2 by duplex 
RNAs that target promoter transcripts near the transcrip- 
tion start site is not a unique observation. Our laboratory 
has previously observed similar results with duplex RNAs 
targeting promoter regions for progesterone receptor 
(20,34,35) and LDL receptor (21). The activating RNAs 
responsible for controlling LDL receptor, progesterone 
receptor and COX-2 expression have little sequence simi- 
larity, and regulate both TATA-less and TATA-box pro- 
moters. For all three gene targets, binding to a nuclear 
promoter transcript was observed. The common 
activating effects by RNAs that are complementary to 
different gene promoters support the conclusion that the 
ability to sequence- specifically recognize the gene target is 
critical for activity. 

Future directions for research 

Our study is a step toward understanding the potential of 
RNAi factors to control gene expression in cell nuclei. 
Much remains to be learned. The presence of AG02 
and other RNAi factors in mammalian cell nuclei con- 
tinues to be controversial despite significant experimental 
evidence that AG02 can be a nuclear factor (26). It is 
essential that experiments be performed to move this 
debate forward. 
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Our work suggests that nuclear RNAi factors can be 
robust regulators of gene expression and that evolutionary 
pressure should lead to other genes being influenced by 
nuclear RNAi. Sequencing of RNAs bound to RNAi 
factors in cell nuclei may help identify likely gene targets 
for endogenous regulation. We have found that validating 
endogenous targets requires a painstaking and resource- 
intensive effort. Because the experiments are demanding, 
using RNA-seq to rank targets is a high priority. 
Experiments that systematically investigate potential en- 
dogenous targets would determine whether nuclear RNAi 
has broad potential for gene regulation and are needed to 
justify the investment required to pursue in vivo 
investigations. 

Our studies are not intended to lay out a complete 
roadmap for how recognition of a promoter transcript 
by a small RNA in complex with RNAi factors can 
affect gene expression. As with protein transcription or 
splicing factors, a full inventory of influential partners 
and elucidation of a step-by-step mechanism will require 
more extensive research and may vary from one gene 
target to the next. Similarly, the importance of histone 
modification, whether it is a cause or consequence of 
gene activation, will be an important goal for future 
research into the molecular details of RNA-mediated 
gene activation. 

We are less optimistic about the potential for using 
vector-based technologies as a tool for investigating 
physiologic systems that, like ours, appear to function in 
cis. For example, we previously observed that 
overexpression of a promoter transcript required for 
RNA-mediated gene activation had no effect on expres- 
sion of the associated mRNA (36). 

The value of vector-based systems for evaluating cyto- 
plasmic miRNA interactions is well-established. In the 
nucleus, however, our data suggest a more complex inter- 
play of promoter transcripts, RNAi factors, histone modi- 
fications and gene loci at the chromatin level. It seems 
unlikely that reporter constructs will often serve as good 
mimics of the endogenous system. Even if effects are 
observed, it is not clear to us that the considerable effort 
needed to set up a reporter system will lead to results that 
will be relevant to the parent gene. At a minimum, the 
disadvantages of moving away from studying the en- 
dogenous gene and the time needed to develop a 
reporter system must be carefully weighed against 
any perceived advantages that a reporter system might 
bring. 

One possible exception would be site-specific mutagen- 
esis of mammalian chromosomes. Such alterations can 
now be achieved with zinc fingers or TALEN technologies 
(49). While beyond the scope of our study, selective mu- 
tagenesis would be important for eliminating any uncer- 
tainty about the direct target for small RNAs that 
modulate transcription. These experiments, however, are 
not trivial. They would require that all chromosomes be 
mutated and that the mutations not disrupt existing tran- 
scription factor binding sites, such as the TATA box site 
adjacent to the miR-589 target at the COX-2 promoter. As 
with the use of vector constructs, there exists a danger that 



altering the endogenous gene might yield results that apply 
only to the mutant construct. 

We have now shown that small RNAs targeted to se- 
quences near transcription start sites can upregulate ex- 
pression of three loci, COX-2 /PLA2G4A, progesterone 
receptor (20,34-37) and LDL receptor (21). Aside from 
any natural regulatory role for small RNAs in mammalian 
nuclei, these data suggest that it might be possible to 
upregulate therapeutic targets. This should be another 
goal for research. 

Our experience suggests that the key step toward thera- 
peutic upregulation is identifying a suitable gene target 
that justifies investing the considerable effort needed to 
find and validate activating RNAs. Once a compelling 
disease target is identified, pursuing experiments similar 
to those outlined in our report should permit a relatively 
rapid yes or no decision about whether RNA-mediated 
activation is feasible and also provide insights into the 
potential mechanism of action. 



CONCLUSION 

Promoter RNA can act in conjunction with small RNAs 
and RNAi factors to organize a multigene network of 
transcriptional regulation and trigger enhanced expres- 
sion. In contrast to miRNAs in the cytoplasm, which 
commonly induce modest changes in gene expression of 
>2-3-fold, we observe a dramatic increase of >25-fold, 
which can be enhanced to > 100-fold when combined 
with other natural activation pathways. 

Our results may be broadly relevant to mammalian gene 
regulation. miRNAs and RNAi factors GW182 and 
AG02 are in the nucleus, providing the macromolecular 
machinery needed for sequence-specific recognition. Non- 
coding RNAs that are associated with chromatin offer a 
recognition point for RNA-mediated control of transcrip- 
tion at specific loci, making the concept a logical starting 
point for understanding the importance of nuclear RNA. 
RNA-mediated gene activation by both fully complemen- 
tary duplex RNAs and partially complementary miR-589 
is robust and potent, making it likely that the pathway 
would be exploited during evolution (26). miR-589 is 
conserved in primates, as are its target sites within the 
COX-2 promoter and the spatial organization of the 
COX-2 and PLA2G4A genes. 

Gene loops are relatively common at gene promoters 
and often extend beyond adjacent genes (50). Together, 
RNA-mediated recognition and ds-acting gene loops 
allow the COX-2 and PLA2G4A loci to form a functional 
operon that can affect the efficiency of eicosanoid biosyn- 
thesis. If RNA-mediated transcription linkages are found 
to be general for other pairs of functionally related 
pathway genes, they would have far reaching implications 
as a regulatory mechanism for sequence- specifically 
controlling the transcription of gene networks. Small 
RNAs that are designed to recognize nascent RNAs can 
achieve potent and selective activation of single genes or 
gene networks and provide a starting point for a new class 
of therapeutic agent. 



10108 Nucleic Acids Research, 2013, Vol. 41, No. 22 



SUPPLEMENTARY DATA 

Supplemental Data is available at NAR Online. 



ACKNOWLEDGEMENTS 

The authors thank Dr Jacob Schwartz, Dr Michael White, 
Dr Elias Corey, Dr Alex Pertsemlidis and Dr Larry 
Marnett for review of the manuscript; Roya Kalantari 
and Dr Liande Li for assistance with experiments; and 
Dr Scott Younger with assistance identifying miR-589. 

FUNDING 

National Institutes of Health (NIH) [GM 73042/DRC, 
GM85080/BAJ]; Welch Foundation [I-1244/DRC]; 
Cancer Prevention and Research Institute of Texas 
[RP120311/BAJ]. Funding for open access charge: NIH. 

Conflict of interest statement. None declared. 



REFERENCES 

1 . The FANTOM consortium and RIKEN Genome Exploration 
Research Group and Genome Science Group. (2005) The 
transcriptional landscape of the mammalian genome. Science, 309, 
1559-1563. 

2. Clark,M.B., Amaral,P.P., Schlesinger,F.J., Dinger, M.E., Taft,R.J., 
RinnJ.L., Ponting,C.P., Stadler,P.F., Morris,K.V., Morillon,A. 

et al. (2011) The reality of pervasive transcription. PLoS Biol., 9, 
el 000625. 

3. Guttman,M. and RinnJ.L. (2012) Modular regulatory principles 
of large non-coding RNAs. Nature, 482, 339-346. 

4. Modarresi,F., Faghihi,M.A., Lopez-Toledano,M.A., Fatemi,R.P., 
Magistri,M., Brothers,S.P., van der Brug,M.P. and Wahlestedt,C. 
(2012) Inhibition of natural antisense transcripts in vivo results in 
gene-specific transcriptional upregulation. Nat. Biotechnol., 30, 
453-459. 

5. Weinberg, M.S. and Morris,K.V. (2013) Long non-coding RNA 
targeting and transcriptional de-repression. Nucleic Acid Ther., 23, 
9-14. 

6. Ogawa,Y., Sun,B.K. and LeeJ.T. (2008) Intersection of the RNA 
interference and X-inactivation pathways. Science, 320, 
1336-1341. 

7. Tsai,M.C, Manor,0., Wan,Y., Mosammaparast,N., WangJ.K., 
Lan,F., Shi,Y., Segal,E. and Chang,H.Y. (2010) Long noncoding 
RNA as modular scaffold of hisone modification complexes. 
Science, 329, 689-693. 

8. Miller,S.B. (2006) Prostaglandins in health and disease: an 
overview. Semin. Arthritis Rheum., 36, 37-49. 

9. Dinarello,C.A. (2010) Anti-inflammatory agents: present and 
future. Cell, 140, 935-950. 

10. Goessling,W., North,T.E., Loewer,S., Lord,A.M., Lee,S., Stoick- 
Cooper,C.L., Weidinger,G., Puder,M., Daley,G.Q., Moon,R.T. 
et al. (2009) Genetic interaction of PGE2 and Wnt signalling 
regulates developmental specification of stem cells and 
regeneration. Cell, 136, 1136-1147. 

11. Lim,H., Paria,B.C, Das,S.K., DinchukJ.E., Langenbach,R., 
TrzaskosJ.M. and Dey,S.K. (1997) Multiple female reproductive 
failures in Cyclooxygenase 2-deficient mice. Cell, 91, 197-208. 

12. TobimD.M., Roca,F.J., Oh,S.F., McFarland,R., Vickery,T.W., 
RayJ.P., Ko,D.C, Zou,Y., Bang,N.D., Chau,T.T.H. et al. (2012) 
Host genotype- specific therapies can optimize the inflammatory 
response to mycobacterial infections. Cell, 148, 434-446. 

13. Wang,D. and DuBois,R.N. (2010) Eicosanoids and cancer. Nat. 
Rev., 10, 181-193. 

14. Vegiopoulos,A., Muller-Decker,K., Strzoda,D., SchmittJ., 
Chichelnitskiy,E., Ostertag,A., Diaz,M.B., RozmanJ., Hrabe de 
Angelis,M., Nusing,R.M. et al. (2010) Cyclooxygenase-2 controls 



enery homeostasis in mice by de novo recruitment of brown 
adipocytes. Science, 328, 1158-1161. 

15. Appleby,S.B., Ristimaki,A., Neilson,K., Narko,K. and Hla,T. 
(1994) Structure of the human cyclo-oxygenase-2 gene. Biochem. 
J., 302, 723-727. 

16. Smale,S.T. (2010) Selective transcription in response to an 
inflammatory stimulus. Cell, 140, 833-844. 

17. Wu,T., Ikezono,T., Angus,C.W. and ShelhamerJ.H. (1994) 
Characterization of the promoter for the human 85 kDa cytosolic 
phospholipase A2 gene. Nucleic Acids Res., 22, 5093-5098. 

18. Ghosh,M., Tucker,D.E., Burchett,S.A. and Leslie,C.C. (2006) 
Properties of the Group IV phospholipase A2 family. Prog. Lipid 
Res., 45, 487-510. 

19. Janowski,B.A., Hu,J. and Corey,D.R. (2006) Silencing gene 
expression by targeting chromosomal DNA with antigene peptide 
nucleic acids and duplex RNAs. Nat. Protoc, 1, 436^43. 

20. Chu,Y., Yue,X., Younger,S.T., Janowski,B.A. and Corey,D.R. 
(2010) Involvement of argonaute proteins in gene silencing and 
activation by RNAs complementary to a noncoding transcript at 
the progesterone receptor promoter. Nucleic Acids Res., 38, 
7736-7748. 

21. Matsui,M., Sakurai,F., Elbashir,S., Foster,D.J., Manoharan,M. 
and Corey,D.R. (2010) Activation of LDL receptor expression of 
small RNAs complementary to a noncoding transcript that 
overlaps the LDLR promoter. Chem. Biol., 17, 1344-1355. 

22. Dodd,D.W., Gagnon,K.T. and Corey,D.R. (2013) Digital 
quantitation of potential therapeutic target RNAs. Nucleic Acids 
Ther., 23, 188-194. 

23. HavelockJ.C, Keller,P., Muleba,N., Mayhew,B.A., Casey,B.M., 
Rainey,W.E. and Word,R.A. (2005) Human myometrial gene 
expression before and during parturition. Biol. Reprod., 12, 
707-719. 

24. Wang,Y., Juranek,S., Li,H., Sheng,G., Tuschl,T. and Patel,D.J. 
(2008) Structure of an argonaute silencing complex with a seed- 
containing guide DNA and target RNA duplex. Nature, 456, 
921-926. 

25. WattsJ.K. and Corey,D.R. (2012) Silencing disease genes in the 
laboratory and the clinic. /. Pathol, 226, 365-379. 

26. Gagnon,K.T. and Corey,D.R. (2012) Argonaute and the nuclear 
RNAs: New pathways for RNA-mediated control of gene 
expression. Nucleic Acid Ther., 22, 3-16. 

27. Robb,G.B., Brown,K.M., KhuranaJ. and Rana,T.M. (2005) 
Specific and potent RNAi in the nucleus of human cells. Nat. 
Struct. Mol. Biol, 12, 133-137. 

28. Liu,J., Hu,J. and Corey,D.R. (2012) Expanding the action of 
duplex RNAs into the nucleus: redirecting alternative splicing. 
Nucleic Acids Res., 40, 1240-1250. 

29. Eulalio,A., Tritschler,F. and Izaurralde,E. (2009) The GW182 
protein family in animal cells: New insights into domains 
required for miRNA-mediated gene silencing. RNA, 15, 
1433-1442. 

30. AzevedoJ., Cooke,R. and Lagrange,T. (2011) Taking RISCs with 
hookers. Curr. Opin. Plant Biol, 14, 594-600. 

31. Lian,S.L., Li,S., Abadal,G.X., Pauley,B.A., Fritzler,M.J. and 
Chan,E.K.L. (2009) The C-terminal half of human Ago2 binds to 
multiple GW-rich regions of GW182 and requires GW182 to 
mediate silencing. RNA, 15, 804-813. 

32. Till,S., Lejeune,E., Thermann,R., Bortfeld,M., Hothorn,M., 
Enderle,D., Heinrich,C, Hentze,M.W. and Ladurner,A.G. (2007) 
A conserved motif in Argonaute-interacting proteins mediates 
functional interactions through the Argonaute PIWI domain. Nat. 
Struct. Mol. Biol, 14, 897-903. 

33. Li,L.C, Okino,S.T., Zhao,H., Pookot,D., Place,R.F., Urakami,S., 
Enokida,H. and Dahiya,R. (2006) Small dsRNAs induce 
transcriptional activation in human cells. Proc. Natl Acad. Sci. 
USA, 103, 17337-17342. 

34. Janowski,B.A., Younger,S.T., Hardy,D.B., Ram,R., Huffman,K.E. 
and Corey, D.R. (2007) Activating gene expression in mammalian 
cells with promoter-targeted duplex RNAs. Nat. Chem. Biol, 3, 
166-173. 

35. SchwartzJ.C, Younger,S.T., Nguyen,N.B., Hardy,D.B., 
Corey,D.R. and Janowski,B.A. (2008) Antisense transcripts are 
targets for small activating RNAs. Nat. Struct. Mol. Biol, 15, 
842-848. 



Nucleic Acids Research, 2013, Vol 41, No. 22 10109 



36. Yue,X., SchwartzJ.C, Younger,S.T., Chu,Y., Gagnon,K., 
Elbashir,S., Janowski,B.A. and Corey,D.R. (2010) Regulation of 
transcription by small RNAs complementary to sequences 
downstream from the 3' termini of genes. Nat. Chem. Biol, 6, 
621-629. 

37. WattsJ.K., Yu,D., Charisse,K., Montaillier,C, Potier,P., 
Manoharan,M. and Corey, D.R. (2010) Effect of chemical 
modifications on modulation of gene expression by duplex 
antigene RNAs (agRNAs) that are complementary to gene 
promoters. Nucleic Acids Res., 38, 5242-5259. 

38. Morris,K.V., Santoso,S., Turner,A.M., Pastori,C. and 
Hawkins,P.G. (2008) Bidirectional transcription directs both 
transcriptional gene activation and suppression in human cells. 
PLoS Genet., 4, el000258. 

39. Lim,L.P., Lau,N.C, Garrett-Engele,P., Grimson,A., SchelterJ.M., 
CastleJ., Bartel,D.P., Linsley,P.S. and Johnson, J.M. (2005) 
Microarray analysis shows that some micro RNAs downregulate 
large numbers of target mRNAs. Nature, 433, 769-773. 

40. Wang,Z., Yong Jeon,H., Rigo,F., Bennett,C.F. and Krainer,A.R. 
(2012) Manipulation of PK-M mutally exclusive alterative splicing 
by antisense oligonucleotides. Open Biol., 2, 120133. 

41. Robbins,M., Judge,A. and MacLachlan,I. (2009) siRNA and 
innate immunity. Oligonucleotides, 19, 89-101. 

42. Liu,J., Carmell,M.A., Rivas,F.V., Marsden,C.G., Thomson, J. M., 
SongJ.J., Hammond,S.M., Joshua-Tor,L. and Hannon,G.J. (2004) 
Argonaute2 is the catalytic engine of mammalian RNAi. Science, 
305, 1437-1441. 



43. Zhang,P., Lee,H., BrunzelleJ.S. and CoutureJ.F. (2012) The 
plasticity of WDR5 peptide-binding cleft enables the binding of 
the SET1 family of histone methyltransferases. Nucleic Acids Res., 
40, 4237-4246. 

44. Wang,K.C, Wang,Y.W., Liu,B., Sanyal,A., Corces- 
Zimmerman,R., Chen,Y., Lajoie,B.R., Protacio,A., Flynn,R.A., 
Gupta,R.A. et al. (2011) A long noncoding RNA maintains 
active chromatin to coordinate homeotic gene expression. Nature, 
412, 120-124. 

45. Li,G., Ruan,X., Auerbach,R.K., Sandhu,K.S., Zheng,M., 
Wang,P., Poh,H.M., Goh,Y., Lim,J., ZhangJ. et al. (2012) 
Extensive promoter-centered chromatin interactions 

provide a topological basis for transcriptional regulation. Cell, 
148, 84-98. 

46. DekkerJ. (2006) The three 'C's of chromosome conformation 
capture: controls, controls, controls. Nat. Methods, 3, 17-21. 

47. Buhler,M., Verdel,A. and Moazed,D. (2006) Tethering RITS to a 
nascent transcript initiates RNAi- and heterochromatin-dependent 
gene silencing. Cell, 125, 873-886. 

48. Cam,H.P., Chen,E.S. and Grewal,S.I.S. (2009) Transcriptional 
scaffolds for heterochromatin assembly. Cell, 136, 610-614. 

49. Gaj,T., Gersbach,C.A. and Barbas,C.F. Ill (2013) ZFN, TALEN, 
and CRISPR/Cas-based methods for genome engineering. Trends 
Biotechnol, 31, 397-405. 

50. Sanyal,A., Lajoie,B.R., Jain,G. and DekkerJ. (2012) The long- 
range interaction landscape of gene promoters. Nature, 489, 
109-113. 



